Numerical Computation of Deflating Subspaces of Skew-Hamiltonian/Hamiltonian Pencils

نویسندگان

  • Peter Benner
  • Ralph Byers
  • Volker Mehrmann
  • Hongguo Xu
چکیده

We discuss the numerical solution of structured generalized eigenvalue problems that arise from linear-quadratic optimal control problems, H∞ optimization, multibody systems, and many other areas of applied mathematics, physics, and chemistry. The classical approach for these problems requires computing invariant and deflating subspaces of matrices and matrix pencils with Hamiltonian and/or skew-Hamiltonian structure. We extend the recently developed methods for Hamiltonian matrices to the general case of skew-Hamiltonian/Hamiltonian pencils. The algorithms circumvent problems with skew-Hamiltonian/Hamiltonian matrix pencils that lack structured Schur forms by embedding them into matrix pencils that always admit a structured Schur form. The rounding error analysis of the resulting algorithms is favorable. For the embedded matrix pencils, the algorithms use structure-preserving unitary matrix computations and are strongly backwards stable, i.e., they compute the exact structured Schur form of a nearby matrix pencil with the same structure.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Computational Experience with Structure-preserving Hamiltonian Solvers in Complex Spaces

Structure-preserving numerical techniques for computation of eigenvalues and stable deflating subspaces of complex skew-Hamiltonian/Hamiltonian matrix pencils, with applications in control systems analysis and design, are presented. The techniques use specialized algorithms to exploit the structure of such matrix pencils: the skew-Hamiltonian/Hamiltonian Schur form decomposition and the periodi...

متن کامل

FORTRAN 77 Subroutines for the Solution of Skew-Hamiltonian/Hamiltonian Eigenproblems - Part II: Implementation and Numerical Results

Skew-Hamiltonian/Hamiltonian matrix pencils λS −H appear in many applications, including linear quadratic optimal control problems, H∞-optimization, certain multi-body systems and many other areas in applied mathematics, physics, and chemistry. In these applications it is necessary to compute certain eigenvalues and/or corresponding deflating subspaces of these matrix pencils. Recently develope...

متن کامل

Robust numerical methods for robust control

We present numerical methods for the solution of the optimal H∞ control problem. In particular, we investigate the iterative part often called the γ-iteration. We derive a method with better robustness in the presence of rounding errors than other existing methods. It remains robust in the presence of rounding errors even as γ approaches its optimal value. For the computation of a suboptimal co...

متن کامل

Numerical Solution of Quadratic Eigenvalue Problems with Structure-Preserving Methods

Numerical methods for the solution of large scale structured quadratic eigenvalue problems are discussed. We describe a new extraction procedure for the computation of eigenvectors and invariant subspaces of skew-Hamiltonian/Hamiltonian pencils using the recently proposed skew-Hamiltonian isotropic implicitly restarted Arnoldi method (SHIRA). As an application we discuss damped gyroscopic syste...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • SIAM J. Matrix Analysis Applications

دوره 24  شماره 

صفحات  -

تاریخ انتشار 2002